Relationship between incommensurability and superconductivity in Peierls distorted 

charge-density-wave systems 
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We study the pairing potential induced by fluctuations around a charge-density wave (CDW) 
with scattering vector Q by means of the Frohlich transformation. For general commensurability 
M, defined as |k + MQ) = |k), we find that the intraband pair scattering within the M subbands 
scales with M whereas the interband pair scattering becomes suppressed with increasing CDW or- 
der parameter. As a consequence superconductivity is suppressed when the Fermi energy is located 
between the subbands as it is usually the case for nesting induced CDWs, but due to the vertex 
renormalization it can be substantially enhanced when the chemical potential is shifted sufficiently 
inside one of the subbands. The model can help to understand the experimentally observed depen- 
dence of the superconducting transition temperature from the stripe phase incommensurability in 
the lanthanum cuprates. 



PACS numbers: 74.20.Fg, 74.25.-q, 74.25.Kc, 74.72.Dn 



I. INTRODUCTION 



It is now well established that some of the cuprate 
high-T c superconductors show a strong depression of 
T c when the doping concentration coincides with an 
integer fraction. This feature was first detected in 
La2- 2: Ba 2: Cu04 |j] where the superconducting transi- 
tion temperature displays a clear minimum near com- 
positions x sw 0.12 w 1/8. Subsequent measurements 
on La2-:rSr K Cu04 (LSCO) also resulted in a detectable 
dip of T c around the same doping concentrations @ al- 
though the reduction is less pronounced in this system 
than in the Ba-doped samples.. Further insight into 
this anomaly was provided by inelastic neutron scattering 
measurements on L^-zSr^CuO^ (LSCO) which re- 
vealed the existence of incommensurate spin correlations 
in this material. More recently it was found that these in- 
commensurate spin fluctuations are pinned in Nd doped 
LSCO oxide compounds ||,|7j leading to the appearance 
of magnetic and charge-order superlattice peaks in the 
neutron diffraction pattern. These measurements there- 
fore provide strong evidence for the occurence of spin- 
and charge-stripe order in the lanthanum compounds 
which becomes static when the lattice undergoes a tran- 
sition to a low temperature tetragonal (LTT) phase as it 
is the case for the Ba- and Nd doped systems. However, 
there is now growing evidence that this structural phase 
transition is not a necessary condition for the occurence 
of static stripe order which recently has also been de- 
tected in LSCO M and oxygen-doped La2Cu04 +!/ ||. It 
should be noted that in case of the LSCO compound, ul- 
trasonic measurements carried out for doping x = 0.12 
show a mode softening below ~ 45 K M], which can be 
understood as a precursor of the LTT structural phase 
transition. This in turn can act as the pinning mechanism 
for the fluctuating stripe order and therefore explain the 
small dip in T c around x « 1/8 for the LSCO material. 



With regard to the origin of incommensurable scat- 
tering it has been concluded that the stripe ordering is 
driven by the charges H|n] since the magnetic peaks ap- 
pear at lower temperatures than the charge-order peaks, 
consistent with an Landau free-energy analysis of cou- 
pled charge- and spin-density order parameters |L2| . For 
completeness we want to mention that incommensurate 
spin correlations have also been detected in the YBCO 
compounds 13 and the existence of charge modulation 
in the Cu02 planes of Bi2212 was proposed early on in 
Ref. @. 

It is therefore naturally to ask wether there may exist a 
connection between the occurence of stripe modulations 
and the superconducting pair formation. A commonly 
held opinion is that stripe correlations naturally lead to 
a depression of T c (as observed in Nd-doped LSCO) due 
to the (nesting induced) reduction of the density of states 
at the Fermi level and the localization of the charge car- 
riers through the pinning potential of a possible struc- 
tural phase transition. However, most interesting results 
have been obtained in the LSCO system where both the 
incommensurate spin fluctuations and superconductivity 
have been studied as a function of doping |D| . It turned 
out that the incommensurability of the magnetic fluctua- 
tions (defined as the deviation of the magnetic peak from 
its antiferromagnetic position) scales linearly with the su- 
perconducting transition temperature T c up to optimal 
doping, i.e. both 8 and T c appear at the same critical 
concentration x» 0.05 and follow the same doping de- 
pendence. This result strongly suggests that incommen- 
surate stripe ordering and superconductivity may arise 
from the same instability as it is proposed in the frame- 
work of a Quantum Critical Point (QCP) scenario Jl6| . 
Within this model the singular scattering induced by the 
critical charge fluctuations would be responsible for both 
the anomalous normal-state properties and the large su- 
perconducting critical temperatures. Locating the QCP 



1 



near optimal doping the phase diagram of the cuprates 
is partitioned in a (nearly) ordered, a quantum critical, 
and a quantum disordered region naturally correspond- 
ing to the under-, optimally, and over-doped regions of 
the phase diagram of the cuprates. 

Presupposing the existence of charge stripe correla- 
tions we want to investigate in the present paper how the 
incommensurability of these modulations may affect the 
pairing potential. For convenience we restrict ourselves 
to the static case using a Frohlich type transformation 
to study the problem how the modification of the elec- 
tron states due to the formation of a charge-density wave 
(CDW) influences the superconducting order parameter 
for given (in)commensurability. It should be mentioned 
that a similar model has already been considered in Ref. 

where the effective electron-electron interaction me- 
diated by phonons has been factorized due to CDW and 
superconducting order parameters. However, this proce- 
dure is not consistent with the formation of a CDW since 
it leads to an unnatural energy cutoff for the particle-hole 
excitations whereas a static CDW involves particle-hole 
scattering processes of all available energies. The correct 
procedure which we will follow below is therefore first 
to diagonalize the CDW mean-field part of the hamilto- 
nian and after that to calculate the effective interaction 
between the new quasiparticles arising from fluctuations 
around the frozen-in CDW. It turns out that this inter- 
action strongly depends on the (in)commensurability M 
of the CDW with the underlying lattice which opens the 
possibility of increasing T c by doping the CDW while 
fixing its periodicity. We want to stress that this mecha- 
nism is not due to the appearance of 1/ Y^-singularities 
in the density of states (DOS) but has its origin in the 
modification of the electron-phonon vertex function. 

In Sec. II we introduce our model hamiltonian and in 
Sec. Ill we will derive the effective quasiparticle interac- 
tion in a two-dimensional CDW system of arbitrary com- 
mensurability using a Frohlich transformation approach. 
Explicit results are presented in Sec. Ill for the simplest 
case of a CDW with modulation Q = (tt,tt). Despite its 
low commensurability it turns out that T c is consider- 
ably enhanced due to the vertex renormalization when 
the chemical potential is located sufficiently inside one of 
the two subbands. We finally summarize our conclusions 
in Sec. IV. 



II. MODEL HAMILTONIAN 

Our investigations are based on the following lattice 
model which describes the coupling of electrons to a dis- 
persionless phonon mode 

H = tkc\ „c Ka + uj b l b i 

k,a q 



+ ^E c U^(fo 9 +6 t _ 9 ) a) 

kq,<T 

where eft destroys (creates) an electron in the state 

k and energy ej,. The operates b q ^ destroy (create) a 
phonon with momentum q. For convenience we take the 
electron-phonon matrix element g and the phonon fre- 
quency ujq to be constant. 

Within this model the transition to a static CDW is 
signaled by a diverging polarizability for a certain wave 
vector Q and at some critical doping concentration. The 
phonon part of the CDW mean-field wave- function is usu- 
ally taken to be a coherent state which means that for 
the selected Q-values the phonon operators are replaced 
by a complex number. 

The coherent state contribution is extracted from the 
hamiltonian Eq. ([!]) by means of a displacement trans- 
formation b q — a q — a q where the operators a q now cor- 
respond to the phonon fluctuations around the coherent 
state \a q ). The CDW order parameter is defined through 
gctq/VN = Aq5 9i q. We have neglected the contribution 
of higher harmonics to the CDW profile which are not 
essential for our present considerations. 

We thus obtain the following hamiltonian 

H = Hcdw + Hpi + Hp (2) 

where 

Hcdw = J2 e * c iL c fc,- + ^I a q| 2 

k,a 9 

~ ( A Q C i+Q,a C k,* + &* Q 4;-Q,* C k,*) (3) 

k,a 

decribes the static CDW scattering with wave vector Q of 
the electrons. The essential term which contains the cou- 
pling of the phonon fluctuations to the electronic degrees 
of freedom corresponds to 

H FL = -r^J2 C l+q,a C kAa q + ol,) (4) 
kq,a 

and 

Hp = ujqY a l a q ( 5 ) 
q 

is the energy of the phonon fluctuations. 

III. EFFECTIVE INTERACTION FOR GENERAL 
COMMENSURABILITY 

For a given k-state in the Brillouin zone we define the 
commensurability M of the CDW as the number of scat- 
tering events in the same direction needed to return to 
the equivalent state, i.e. |k + MQ) = |k). It is clear 
that the larger M the higher is the 'incommensurability' 
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of the CDW with the underlying lattice. In order to keep 
the comprehensibility of our presentation we consider in 
the following a two-dimensional system where the CDW 
scattering vector Q = (Q x ,0) is along the x-direction. 
Within the derivation presented below other orientations 
would only change the k-state labeling. The reduced Bril- 
louin zone (BZ) then is defined by —tt/M < k x < tt/M 
and — 7r/a < k y < tt / a and the electron operators in the 
reduced zone read as 



C n (k,a) = C fcx _ sgn (fc x )(n--l)Q x ,fc B ,<T 



(0) 



where the same definition holds also for the phonon op- 
erators. The signum function has been introduced in 
order to obtain symmetric subbands in the reduced zone 
with respect to k x . Performing a unitary transformation 
c n (k, er) = X)m=i A nm f m (k, a) the mean-field part of the 
hamiltonian Eq. (J3J) can be diagonalized 

Hcdw = E E E n(k)fl{k, a)f n (k, a) (7) 

k,a n 

where E n (k) denotes the dispersion of the n-th sub- 
band in the reduced BZ. It can be further shown that 
the condition E n (k x ,k y ) = E n (—k x ,k y ) requires for 
the transformation matrices the property A nm (k x , k y ) = 
A nm (~kx,k y ) = A mn {— k x , ky)- 

The coupling term Eq. (|J) rewritten in terms of the 
new quasiparticle operators reads as 

Hfl = -j=^2^2r stm (k,q)f s (k + q,a)f t (k,a) 

kq,cr stm 

x (a m (q) +al n (-q)) 
where we have introduced the vertex function 



(8) 



Alr, s {k + q)A n>t {k) 



and the index TV = Af(n, m, k, q) is defined below. Un- 
fortunately at this point the notation is getting complex 
both due to the addition of momenta in the signum func- 
tion and the scattering of electrons between neighboring 
reduced zones. For this reason the following definitions 
have been made in Eqs. (|[ [)]) 

For sgn(fc :E ) = sgn(q x ) and \k x + q x \ < tt/M : 
k + q = k + q; TV = n + m - 1 (10) 

For sgn(fc K ) ^ sgn{q x ) and \k x + q x \ < tt/M : 
k + q = k + q; TV = sgn(|fc x | - \q x \)(n - m) + 1 (11) 

For sgn(k x ) = sgn(q x ) and \k x + q x \ > tt/M : 
k + q = k + q — sgn(k x + q x )Q x \ TV = n + m — 2 (12) 

The coupling term Eq. (JsJ) can now be eliminated in 
first order by means of a standard Frohlich transforma- 
tion (see e.g. Q), i.e. we expand the hamiltonian Eq. 



(g) in a commutator series up to second order in the ver- 
tex F 



H = e~ s He s = H+[H,S} + ]-[H, S],S] 
with S being specified by the condition: 
[Hcdw + Hp, S] = —Hpi 
Representing S in the form 

S 



(13) 



(14) 



1 



EE ft 



x \P- 



fl(k + q,a)f v (k,a) 

kq,cr uvw 

(k,q)a w (q) + Pu VW { k ^) a w{-Q)] 



the requirement Eq. (^) is satisfied by 

r uvw {k, (?) 



E u {k + q) -E v (k)±L>o 



(15) 



(16) 



We thus finally obtain a transformed hamiltonian 
which, among terms describing the renormalization of 
phonon and quasiparticle energies, contains an effective 
quasiparticle interaction. The relevant contribution for 
superconductivity consists of those terms which describe 
the scattering of a zero momentum pair between the sub- 
bands of the CDW system and reads as 



kq,(T stm 



[E s {k + q)-E t {k)Y-ul 



x fHk + q,a)ft(-k - q,-a)f t (-k,-a)f t (k,a) (17) 

Obviously this expression reduces to the well-known 
phonon mediated electron-electron interaction in the 
limit An — » when the transformation matrices 



(9) are given by A nm (k) 



# n ,m an d the vertex con- 
tribution becomes independent of the band indices 
J2 m \^stm(k, q)\ 2 — > 1. However, in case of a well de- 
velopped CDW the transformation takes the form 



A np (k) 



1 



P 



exp (inpsgn(k x )) 
M 

tt/M ...tt/M 



(18) 



which for the vertex function leads to T s t m {k,q) = 
exp(icj)(m))6s,t i-e. the phonon branch-index m only en- 
ters through a phase <fr(m) depending on the three ranges 
defined in Eqs. ( |Io| - |l2] ). As a consequence the vertex con- 
tribution 



(19) 



now restricts the pair-pair scattering to the individual 
subbands and the effective interaction becomes propor- 
tional to the commensurability M of the CDW. Since 
the number of states within each subband decreases with 



3 



the commensurability as 2N/M it is obvious that the 
vertex renormalization can only lead to an effective en- 
hancement of the quasiparticle interaction as long as the 
cutoff range of its attractive contribution lies within one 
of the subbands. For the present two-dimensional model 
with a phonon mediated interaction this means that for 
M < B/u>d the CDW formation can lead to a substan- 
tial enhancement of the superconducting order parameter 
where B is the bandwidth and lod denotes the Debye fre- 
quency as a characteristic energy scale for the cutoff. Of 
course the considered limit of a large CDW order param- 
eter is rather academic since then the quasiparticles be- 
come localized and superconductivity will be suppressed. 
However, the essential point is that with the onset of 
CDW scattering the intraband pairing is amplified with 
a simultaneous decrease of the interband pairing ampli- 
tude. Both processes do not compensate since most of 
the interband contribution to the interaction falls out of 
the range of attraction. 









64 ; 




32 


s 



0.1 0.2 

|A| 

FIG. 1. Intraband vertex contribution to the effective inter- 
action M a , a = J2 m |F, sm (fc,q)| 2 for Q = (37r/4,0), (5tt/8,0), 
(97r/16,0), and (197r/32) corresponding to the commensu- 
rabilities M = 8, 16, 32, and 64 as indicated in the figure, 
k = (tt/80, it), q = (-vr/40, 0) 

In the remainder of this section we consider as an ex- 
ample the CDW induced vertex renormalization for holes 
with energy dispersion 

efc = 2ti (cos k x + cos k y ) + 4^2 cos k x cos k y 
+ 2*3 (cos 2k x + cos 2k y ) 

where we use the hopping integrals t\ — 45meV, *2 = 
5meV, 7j 3 = 16meV [^0|. For small and moderate hole 
concentrations this dispersion gives rise to an open (i.e. 
centered around (tt, tt)) Fermi surface as observed in var- 
ious photoemission experiments in the high-T c cuprates 
(see e.g. [|l9| and references therein). Moreover one finds 
that in the 'underdoped' regime this model leads to a 
polarizability which is peaked along the k x - and k y - axis 



of the Brillouin zone repectively. As a consequence the 
system becomes unstable towards a vertically oriented 
CDW when the electron-phonon coupling constant ex- 
ceeds some critical value. In the following we consider 
a scattering vector Q ss (0.6-7T, 0) which in our model is 
realized for hole concentrations S as 0.25. 

Fig. 1 displays the intraband vertex contribution to 
the effective interaction M SiS (k,q) = J2 m \^ssm(k, q)\ 2 
for different commcnsurabilities, i.e. the scattering vec- 
tor Q as (0.67T, 0) is approximated by increasing fractions 
(n/m, 0) as indicated in the figure capture. For each case 
we have selected a subband which covers the Fermi level 
and the phase of the CDW order parameter is fixed to 
4> = 7r/3. The general behavior of M StS (k,q) is consis- 
tent with the limiting values discussed above. Starting 
from M S)S (fc, q) = 1 one observes a strong increase for 
some critical value of the CDW order parameter |A| cr it 
and a saturation to the maximum commensurability in 
the limit of large |A|. The value of |A| cr j( strongly de- 
pends on the selected band and especially on the values 
of k and q. When these states are choosen such that 
their scattered states (i.e. k + Q, k + 2Q ... and q + Q, 
q-l- 2Q ...) do not differ much in energy, the onset of the 
increase of M s>s (k, q) occurs for smaller | A| as in case of 
large energy differences between the scattered states. 




0.05 0.1 0.05 0.1 



|A| |A| 

FIG. 2. Vertex contribution to the effective interaction 
M s ,t = ^2 m \T s tm(k,q)\ 2 for two energetically neighbored 
subbands of a CDW with scattering vector Q = (57r/8,0) 
(commensurability M=16). The solid lines correspond to the 
intraband scattering and A/2,2) whereas the dashed line 

refers to the scattering between the subbands (Mi^). a) and 
b) are plots for two different phases (j> = 0, 7r/7 of the CDW 
order parameter A = |A| exp(i^). 

Another point not discussed so far is the influence of 
the phase of the CDW order parameter. In Fig. 2 we 
show the vertex contribution M Stt for two energetically 
neighbored subbands of a CDW with scattering vector 
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Q = (57r/8,0) (commensurability M=16) and two differ- 
ent phases <fi — 0, ir/7. In case of <f> — not only the 
intraband scattering is enhanced with increasing A but 
also the scattering between the bands due to their large 
degeneracy. As a result the curves for the intraband con- 
tribution approach half the value of the commensurabil- 
ity only in the limit of large |A|, in favor of a similar 
behavior of the interband renormalization. As shown 
in Fig. 2b this interband scattering can be suppressed 
by choosing a different phase which essentially restores 
the 'orthogonality' between the quasi-degenerate states. 
However, for the quasiparticle pairing the phase of the 
CDW order parameter is not relevant since it essentially 
redistributes the pair-scattering between different chan- 
nels each having the same pairing amplitude. 



IV. CDW WITH COMMENSURABILITY M=2 

We will now discuss the consequences of the renormal- 
ized vertex to superconductivity in more detail studying 
the simplest available CDW system in two dimensions 
with scattering vector Q = (7r, it) and electron dispersion 
(k = —£k+Q = —B/2(cosk x + cosfc y ), corresponding to 
nearest neighbor hopping of particles on a square lattice. 
The hamiltonian Eq. (|^) then can be easily diagonalized 
and one obtains 

#cdw = V EkifliK a)h(k, a) ~ &{k, a)fc(k, a)] 



NuJ A 2 



(20) 



where the dispersion of the CDW system is given by 
E k = v/A£ DW + i{. 

Following the procedure described above we then de- 
rive the pairing hamiltonian 



wo 



kq,a i,j 



[E(k + q)-E(k)] 2 -uJ 2 



—V for k-states out to a cutoff energy ujd- The sym- 
metry of the SC gap is then completely determined by 
the symmetry of T? ■ (kq) and takes the form 



(fc) = A + a 



A 2 
Ek 



(23) 



0.3 



* 0.2 
E 



0.1 




x fl(- k -q,-a)f}(k + q,a)f j (k > a)f j (-k,-a) (21) 
where the vertex contribution reads 



i± 



'CDW 

EkEk+q 



(22) 



and the +(— ) sign corresponds to i = j(i ^ j). The pair 
scattering between valence and conduction band thus is 
decreased with increasing CDW order parameter in favor 
of an enhancement of the intraband pairing amplitude. 

We have calculated selfconsistently the CDW- and 
SC order parameters within a BCS decoupling scheme 
and furthermore have used a BCS-type approxima- 
tion for the interaction in Eq. (pl|), i.e. replacing 
g 2 u)o/[E(k + q) — E(k)] 2 — by a constant potential 



0.8 0.85 0.9 0.95 
concentration 

FIG. 3. CDW- and SC order parameters as a function of 
doping for a Q = (tt, n) CDW system. Dashed line: CDW 
gap (Acdw), Dotted line: constant part of the SC gap (Ao), 
Solid: Maximum value of the total SC gap (A™^(k)). For 
convenience the calculation has been performed for a con- 
stant DOS p(lj) = l/(2B)e(B 2 -lu 2 ). Parameters: B=0.8eV, 
V=0.25eV, lud =0.04eV. 

Fig. 3 displays the results as a function of doping (con- 
centration n—\ corresponds to half-filling) for a selected 
parameter set. The onset of CDW scattering leads to a 
l/v^o; contribution to the DOS around the CDW gap- 
edges. As a consequence also for concentrations n < 1 
the DOS at the Fermi energy is increased leading to an 
enhancement of the constant part of the SC gap. Be- 
sides this well known effect it can be seen that the ver- 
tex renormalization induces a considerable increase of the 
SC order parameter due to the amplification of the in- 
traband pairing. Additionally superconductivity is sup- 
pressed when the concentration leads to perfect nesting 
of the Fermi surface (i.e. n = I). In the small coupling 
limit one can obtain an explicit expression for the vertex 
contribution to the SC gap under the assumption that 
the chemical potential [i is sufficiently inside the lower 
subband (/1 < — Acdw -wb). We find 



V 



A 2 



, rA ln[ 
4B\n\ 1 4 



(IA4 + 1/ M 2 )] 



(24) 



where B denotes the electronic bandwidth. Note that 
a > since the expression has been obtained under the 
assumption An < ujd, ImI- 
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V. CONCLUSION 

To summarize, we have shown that the vertex renor- 
malization can lead to a substantial enhancement of the 
superconducting order parameter in a Peierls distorted 
CDW system when the chemical potential is located suf- 
ficiently inside one of the subbands. Since for a nest- 
ing induced instability the Fermi level always falls in the 
DOS minimum between two subbands it is clear that this 
mechanism can only work when the CDW modulation is 
fixed by some external (e.g. lattice-) potential while the 
carrier concentration is varied. This in turn, however, is 
not in favor of superconductivity due to the localization 
of the charge carriers. The more promising alternative 
would be a nesting independent CDW instablity as re- 
alized for example in a frustrated phase separation sce- 
nario In this case the CDW formation occurs due 
to the interplay between long-range Coulomb forces and 
a phase separation instability and is not related to any 
Fermi wave vector hp- 

Our considerations provide a possible explanation for 
the linear relationship between T c and incommensurabil- 
ity as observed in the underdoped lanthanum cuprates 
JL5|. In these systems the incommensurability M is re- 
lated to the deviation of the magnetic scattering vector 
Q = (tt/M,tt) from (tt, 7t), i.e. upon doping towards 
rpmax Qne observes an increase of the incommensurabil- 
ity M. According to our findings this in turn induces an 
enhancement of the effective interaction V e f / ~ M which 
causes T c to scale with the incommensurability M. Of 
course we have used a very simplified model restricting 
ourselves to static CDW formation and neglecting spin 
degrees of freedom. However the main outcome of our 
investigations, namely the scaling of the effective interac- 
tion with the commensurability, should hold also in case 
of dynamic CDW scattering. Indeed there is experimen- 
tal evidence that in the La2- a; Sr : i;Cu04 system, where no 
pinning is provided by the LTO-LTT transition, stripes 
fluctuate very slowly plj allowing for an order parameter 
description in principle also in this case. 
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